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Abstract 

A method based on wavelet techniques has been developed and applied to back- 
ground rejection in the data of the IGEX dark matter experiment. The method 
is presented and described in some detail to show how it efficiently rejects events 
coming from noise and microphonism through a mathematical inspection of their 
recorded pulse shape. The result of the application of the method to the last data 
of IGEX is presented. 

PACs: 95.35+d: 14.80.Mz 

Keywords: Dark Matter; Low background; Underground Detectors; Wavelets; 
Noise Rejection 

1 Introduction 

Experiments searching for rare event phenomena, like those looking for non-baryonic 
dark matter particles (axions or WIMPs), supposedly filling a substantial part of 
the galactic haloes, or those looking for the neutrinoless double beta decay have 
become of most relevance in Particle Physics and Cosmology. To find the nature 
of Dark Matter, the essential component of the Universe, which in its non-baryonic 
form might be (according to recent developments) responsible of ~25-30 % of the 
dark energy/matter budget of a flat universe, is one of the big challenges of modern 
Cosmology. On the other hand, the old subject of nuclear double beta decay which 
has been the objective of many experimental efforts as an unique tool to explore the 
nature of the electron neutrino and to verify the lepton number conservation, has 
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become recently of even greater importance after the confirmation that the neutrino 
has a non-zero mass. In fact, the existence of neutrinoless double beta decay may 
provide information on the absolute scale mass of the neutrino, its mass pattern and 
possibly on CP violation in the leptonic sector if the sensitivity of the experiments 
is able to reach the level of a few meV for the Majorana neutrino effective mass 
upper bound. 

The signals to be expected in the two selected examples of rare event physics, 
have in common their low rates and so the main challenge to achieve the required 
sensitivity is to disentangle the signal from the various sources of background, which 
largely hide the expected signals. Consequently, the first basic requirement of these 
types of searches is to provide an ultra-low background environment, as well as to 
design specific techniques of shielding and background rejection. 

In the particular case of dark matter searches, the expected signal is concen- 
trated in the lowest part of the spectrum fl] often just above the threshold of the 
experiment. When the effective threshold of these experiments is fixed by the elec- 
tronic baseline, random fluctuations of this baseline -electronic noise- can populate 
the energy range of interest with noise events. Actually this electronic noise becomes 
a primary source of background which is of special concern. 

In this paper we describe a method which has been developed and successfully 
applied to the last data set of the IGEX dark matter experiment. In the section 
121 we will describe briefly the IGEX experiment, just to understand the needs that 
led to the method presented in this paper. The details of the experiment can be 
found in ref. [21 ^ which contain also the scientific results of IGEX. In sections 3 
and 4 the wavelet method itself and its calibration are described. We finish with 
the conclusions in section 5. 

2 The Experiment 

The IGEX experiment |21I5) optimized for detecting ''^Ge double-beta decay, has 
been described in detail elsewhere. One of the IGEX detectors of 2.2 kg, enriched 
up to 86 % in ^^Ge, is being used to look for WIMPs interacting coherently with the 
germanium nuclei. Its active mass is ~ 2.0 kg, measured with a collimated source of 
^^■^Eu. The full-width at half-maximum (FWHM) energy resolution is 2.37 keV at 
the 1333 keV line of ^'^Co, and the low energy long-term energy resolution (FWHM) 
is less than 1 keV at the 46.5 line of the ^^°Pb. 

We refer to papers jS] and [31 where the latest results of the experiment (regard- 
ing dark matter searches) are presented and where the aspects related to the set-up 
of the experiment, shieldings, the data acquisition system, etc. are described. The 
threshold of the experiment is about 4 keV and the raw background registered in the 
region just above threshold is of a few tenths of counts per kg and day and per keV. 
The background and threshold achievements provide the best WIMP-nucleon cross- 
section versus WIMP mass (cr,m) exclusion plot for spin- independent interactions 
ever obtained by an experiment with no mechanism of nuclear/electron recoils dis- 
crimination for background rejection. As stressed in jHllS], the background level at 
low energies is crucial to set the sensitivity of the experiment and this is illustrated 



2 



by the fact that a moderate reduction of the low energy background would lead 
to a substantial improvement in the final exclusion. To deal with the worrysome 
noise component in this energy region trying its minimization, we developed the 
analysis method described in the present paper. The acquisition system for low 
energy data is implemented by splitting the preamplifier output pulses and routing 
them through two Canberra 2020 amplifiers having different shaping times. These 
amplifier outputs are converted using 200 MHz Wilkinson-type Canberra analog- 
to-digital converters (ADC), controled by a PC through parallel interfaces. The 
energy spectrum is built using the ADC output of one of the lines. In addition, 
the output pulse of one of the amplifiers is digitized and recorded for each event 
using a 800 MHz LeCroy 9362 digital scope. The digitized pulse encompasses a time 
window of 500 /xs with a sampling resolution of 1 fis. This recorded pulse shape is 
analyzed off-line with the method described below and is formally referred as f{t) 
in the following. 



Wavelets are mathematical tools that have been proven to be extremely useful in a 
number of different types of analysis ,6,. In particular, they have been applied to 
the recognition and extraction of noise from certain sets of data. From a mathe- 
matical point of view, the so-called wavelet transforms can be viewed as extensions 
of the Fourier transforms, but with the addition of sensitivity to local features of 
the studied function. In fact, the set of functions which play the role of the e*'^* in 
the Fourier transform are now 'i/'a,fe(^^) derived from an original localized null- area 
wavelet function 'tp{t) by means of translations and contractions parameterized re- 
spectively by b and a. Therefore, for a fixed wavelet function ip{t), the biparametric 
wavelet transform [Wjpf]{a,b) of a given function f{t) can be defined as: 



The above mentioned condition of null-area for ip{t), or admisibility condition, is 
derived from the requirement that no information loss must occur in the transforma- 
tion (i.e. it must be invertible). Apart from this requirement, other properties are 
desirable in certain contexts, like continuity, gaussian shape, etc. Without entering 
into the mathematical details of the wavelet theory [H], and following a practical 
approach, the "mexican hat" wavelet function ip{x) was chosen because it fitted 
most suitably our requirements: 



As stated before, f{t) is the pulse shape for each event as recorded by the 
data acquisition system. Following expression the wavelet transform function 
[W^f]{a,b) can be calculated numerically for each event. The resulting function 
is continuous due to the properties of the wavelet function used. So the relative 
maxima (aj,6j)j (where i goes from 1 to N, being N the total number of maxima) 



3 The method 
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can be calculated by means of an appropriate numerical algorithm. Let's order the 
maxima (oj, bi) according to the magnitude of the wavelet transform in that maxima, 
iOi = [W^f]{ai, bi), such that the first maximum wi is the highest one. In some sense, 
these maxima contain the information of the pulse. For instance, the position of the 
highest one in the (a, b) plane is determined by the width and position in the time 
window of the event pulse itself, if present, and the others (oj , bi)i=2...N are the result 
of random fluctuations of the baseline along the time window encompassed by f{t). 
That is illustrated in figure ^ where in the lower part a sample event pulse is shown 
(to give an idea of the vertical scale, the pulse correspond to an event of about 4.5 
keV, i.e. just above our threshold). In the upper part of the figure the 2-dimensional 
wavelet transform of the same pulse is presented, lighter regions corresponding to 
higher values of the \W^f]{a,b). The relative maxima are also marked with darker 
spots, the first one being clearly identified and corresponding to the true pulse, and 
the others corresponding to random fluctuations of the baseline. 

The key idea of the method is based on the comparison between the magnitude 
of the first maximum cji with the distribution of the others {u)i}i=2...N- In fact, 
the mentioned distribution has been checked to be well modelled by a decaying 
exponential (this fact being related to the kind of wavelet function we are using). 
Figure 121 shows such a distribution for one event. 

Now, the method consists of constructing the exponential distribution D{w) = 
j^^-aw -^j^jg}^ follows the baseline induced maxima in every event pulse shape and 
evaluate the parameter Pi of the first maximum, defined as the normalized area of 
the distribution D{w) which lies above the first maximum wi: 



This parameter Pi is the final outcome of the method for each event, and gives 
us, to a certain extent, an indication of the probability of the main pulse to have 
been generated randomly by the fiuctuations of the baseline, being smaller for pulses 
farther from a noise-like pulse. Therefore, if Pi is large enough, the event must be 
considered a noise fiuctuation and should be rejected, otherwise, the event is kept. 
The reference value to decide if Pi is large enough or not must be fixed by calibration 
of the method, and is explained in the following section. 



In order to fix the rejection criterion this method was applied to a calibration set 
of data. The calibration was performed with a ^^Na source which was introduced 
inside the shielding through a teflon probe producing a rate of events much higher 
than the standard background. The source provides an almost pure sample of 
"true" events, i.e. events not coming from noise. The result of the method applied 
to this set of data is illustrated in figure 01 where the parameter Pi (obtained as 
previously explained) versus energy is presented for each event of the calibration 
set. As expected, this parameter shows a clear correlation with the event energy. It 
is interesting to see that the same plot for the parameter P2, defined similarly as Pi 



Pi 




= e 



(3) 



4 Calibration 
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but for the second maximum W2, shown in figure^lacks completely this correlation, 
as this second maximum is related to some fluctuation of the baseline other than 
the true pulse. The same type of plot obtained for a background set of data is 
shown in figure [3 In this plot two different populations are clearly visible. One of 
them follows the same energy dependence than that of the calibration data and the 
other does not. This second type of events are identified as noisy fluctuations of 
the baseline that have triggered the acquisition system but no relevant pulse above 
threshold is present. From these plots a criterion of Pi < 0.01 can be defined to 
distinguish the two populations of noise and data. Due to the limited statistics of 
the set of data available, it is hard to quantify the energy dependence of the loss 
of efficiency of the cut. However, from the presented plots one can state that it is 
negligible for events which lay a little above our threshold, whereas only in the very 
first keV above threshold it could be only a few percent. 

Moreover, this method provides another way of rejection. The calculation of the 
first maximum of the wavelet transform provides the width of the corresponding 
pulse: ai. This parameter is tightly constrained for true pulses due to the shaping. 
As can be seen in the figure El for our present set-up all of the low energy pulse 
widths are between 7 and 16 fis. This fact can be used to reject events with pulse 
widths outside this range. This argument allows to reject a small population of 
events which were not rejected by the previous method, fact that is illustrated 
in figure 13 They correspond to large lower-frequency fluctuations of the baseline 
presumably connected with microphonism. Due to the fact that they are not strictly 
random fluctuations in the sense previously considered, they could not be rejected 
before. The loss of efficiency of this cut is completely negligible. 

5 Conclusions 

We have presented a method based on wavelet techniques which is able to efhciently 
reject electronic noise and microphonism from the low energy -just above threshold- 
region. The method analyzes the recorded pulse shape of the event to extract 
information about the pulse height in comparison to the fluctuations of the baseline 
to set if the pulse can be assumed to be randomly generated by the fluctuations. 
This method was succesfully applied in the last data of the IGEX dark matter 
direct search. In figure |H1 we show the low energy spectrum obtained with the last 
published 80 kg- day of data of IGEX showing the effective rejection achieved 
with the present method. About half of the background between 4 and 7 keV and 
about 20% between 7 and 10 keV was identified as noise and eventually rejected. 
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Figure 1: In the lower part a sample pulse corresponding to an event of about 4.5 keV 
is shown. In the upper part the 2-dimensional wavelet transform of the same pulse is 
presented. Lighter regions correspond to highest values of the transform function. The 
relative maxima are marked with dark spots. The highest maximum is clearly related 
with the presence of the true pulse, the others being caused by random fluctuations of 
the baseline. 
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Figure 2: Distribution of the maxima a;, {i going from 2 to A^) of the wavelet transform 
for a sample event and their fit to an exponential. 
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Figure 3: Scatter plot showing the parameter Pi assigned to each event by the wavelet 
technique (described in the text) versus energy for a calibration set of data. 
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Figure 4: Scatter plot showing the parameter P2 related to the second maximum of the 
wavelet transform versus the energy of the event (for a calibration set of data) . 
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Figure 5: Same as figure El but for background data. The populations of noise and data 
are well separated above 4 keV. 
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Figure 6: Distribution of pulse widths for a calibration set of data considering only events 
from 4 to 40 keV. All of them are in the range 7 to 16 //s. 
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Figure 7: Distribution of pulse widths for a background set of data (considering only 
events from 4 to 40 keV), before and after the rejection based on the criterion on the 
parameter Pi. It can be appreciated some events with large width (above 16 fis) which 
are not rejected by such criterion. 
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Figure 8: Low energy spectrum of the last published 80 kg- day of data of IGEX before 
(thin line) and after (thick line) the application of the wavelet rejection method. 
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